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EKXL: A DYNAMIC POISSON-SOLVER FOR SPHERIC ALLY-CONVERGENT 
INERTIAL— ELECTROSTATIC CONFINEMENT SYSTEMS 


I. HISTORY AND BACKGROUND OF THE EKXL CODE 

The EKXL code, and each of its predecessors, is a 1-D radially-dependent Poisson-solver. 
It gives static solutions for potential and density distributions, given edge density 
specifications for these parameters. It is a "point” model in that spatial variation of 
externally— imposed fields is possible, but spatial dependence of inter-particle interaction 
phenomena can not be accommodated as an integral part of code operations. 

The original code, IACCEL2/EACCEL2, was developed by Ensley 1 * 2 and provided closed 
analytic solutions in integral form using Struve and Bessel functions. This code employed 
Gaussian distributions for transverse energy and momentum. Solutions were limited, 
however, to electron currents of 1-10 A and core densities of Ixl0 9 -lxl0 10 1/cm 3 , as 
indicated in Figure 1. 

The first modification to Ensley’s work was ACCEL, developed by Mission Research 
Corporation 3 (MRC). This version of the code utilized a more efficient algorithm for 
converging a solution to Poisson’s equation, resulting in significant run-time savings. New 
physics elements were also added in order to include tangential current densities where 
previously only radial current was used. As shown in Figure 1, with the ACCEL 
improvements solutions were achieved with electron currents of 100-500 A and core 
densities of lxlO 11 — lxlO 12 1/cm 3 . 

What followed was a new code, XL, also developed by MRC 4 , that was utilized to its 
capacity by EMC 2 to study the Polywell tm /HEPS system®. This code employed a 
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differential numerical analysis procedure and matrix inversion solution method. The 
distribution functions for transverse energy and momentum were changed to square 
functions to simplify the calculation. Features added at this time were adjustable precision 
and an adaptive grid. This code achieved much higher electron recirculating currents, 
lxlO 9 A, and core densities, Ixl0 13 -lxl0 14 1/cm 3 . However, the core density barrier at 
lxlO 14 1/cm 3 could not be surpassed, making any studies of Polywell tm reactor densities 
impossible. Work was then underway at MRC on the next revision, KXL. 

H. THE ADIABATIC KXL CODE 

KXL uses the same mathematical engine developed in the XL code and a series of 
intermediate "steady-state" solutions over time to reach the final steady-state condition 
for the given input parameters. Selected parameters, such as input current and the length 
of the time step, can be adjusted as the calculation proceeds. As indicated in Figure 1, this 
method and some of the improvements discussed below allow the code to reach much 
higher core densities than were possible with XL. 

The "trapping" formula, F trap , allows feedback on the current density to adjust the input 
current. The algorithm adjusts the current based on the core potential to prevent the 
anode from maximizing, due to too many ions in the core, causing the calculation to be 
unsuccessful. A target anode potential and maximum beam current are input and the code 
adjusts the beam current up to tflft maximum in an attempt to obtain the target anode 
height. The input variable is V 0 and it represents the target well depth. V 0 can be 
different for each species, allowing separate control over each beam. For example, a V Q for 
the electron beam greater than the maximum well depth would set F trap equal to one and 
allow the electron beam to remain on at maximum current for the entire run while the code 
adjusted the ion current wi'hin the indicated range. The anode control formula is: 
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where a = 


( 1 ) 


1+e 


±.03 <p 


determined by the charge on the beam. 


If the maximum available ion current is too high, however, the F trap formula cannot turn 
the current down enough. In effect, there is gun "leakage" and the ions build up (due to 
their recirculation) and overwhelm the core. When using a large ion current, it still needs 
to be within a reasonable range of the final gun current for the code to run successfully. 


The adaptive grid allows for a specified percentage of grid points to be redistributed by the 
code, as needed, to resolve high density and/or high gradient regions. A radial weighting 
factor ("R" weighting) was first used to ensure an adequate number of grid points near the 
outer radius, a region where the density changes rapidly. However, this weighting of the 
outer region did not allow sufficient grid points in the core to resolve high density cases. 
The adaptive gridding was then changed to weight on density as opposed to radius. Once 
this was done, the cases began to run more consistently and the core resolution was greatly 
improved. Figure 2a shows a KXL run with the "R" weighting and Figure 2b shows the 
same case with the "n" weighting. Note the improved resolution of the core region in 
Figure 2b. 

The KXL code was run on the Cyg computer to reach the high density solutions that 
appeared unattainable on the PC, in an attempt to fill in high density points in a study of 
core density vs. electron current. Examining this data, however, it became apparent that 
the core densities were much lower than expected from earlier calculations. Numerous 
variations cf input variables were tried in an attempt to increase the densities. In the 
mean time, an explanation was sought as to why the densities were lower than predicted, 
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and attempts were made to model this behavior. Also, the slope of core density vs. 
electron current was expected to be linear for constant electron confinement. As Figure 3 
indicates, it was not linear until the electron current reaches 1000 A. 

The explanation for the change in slope at 1000 A was found, but without illuminating the 
too-low density problem. The electron density is increasing linearly the entire time, as 
expected. It is only the ion density that changes, due to the anode height. The £n required 
to form an anode is approximately 6n = (n. - n^) = 3.32xl0 9 (E anode /r^), for n in 1/cm 3 , 

E in keV and r Q in cm. The electron currents below the linear slope corresponded to 
densities where £n was a significant portion of the core density. As core density moves 
above this density, 6n becomes less significant, the effect disappears, and the slope finally 
becomes linear. By plotting this data for varying anode heights it is apparent that the 
curve approaches the linear condition as the anode height approaches zero (Figure 4). 


In an effort to understand the too-low core densities, attention turned to the transit times. 
The code was modified to print the transit time at each step and this information was 
plotted, revealing that the code-calculated transit times were much too small to 
correspond to reality (Figure 5). Several changes were suggested and tried to improve the 
transit times, however it was then found that altering the transit time calculation had no 
effect on the final solution. This led to a suspicion that the dn/dt equation, where n is the 
edge density, that was used to link the intermediate steady-state ("snapshot") f jlutions 
might be incorrect. In the KXCTode this equation (for electrons) had the form: 


dn 
e 

dt 


( 2 ) 
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which gives the final steady-state relation n. n .Gj = The basic premise of this 
time-dependent coupling equation (2) is, in fact, false. The correct formulation is one that 
connects time-dependence of total electron populations in the system between successive 
static solutions from the KXL code. This equation is: 


dn 47rR 2 n. . v. . n 

„ inj inj 

dt 4/3 ttR 3 t c 


( 3 ) 


However, since the codes are all written to solve only one-dimensional boundary-value 
problems, it was necessary to recast this equation (3) in terms of edge (boundary) densities. 
And, other physics phenomena of importance to system operation also needed to be added 
to improve the modeling capability of the code. 

Taking all of these factors together, EMC 2 then began development of the EKXL code, 
using the mathematical engine of KXL but now including all of the appropriate 
ion/electron/plasma physics analytic forms appropriate to and required for 
spherically-convergent inertial-electrostatic confinement systems. 


m. THE EKXL CODE 


Change in Edge Density with Respect to Time 

The first step in customizing EKXL to the Polywell tm system was to derive a new, correct 
dn/dt equation. This equation is the key to the linkage of intermediate steady-state 
solutions over time, and numerical tests have shown solutions to be very sensitive to errors 
in this equation. The change in the total number of particles in the system is simply the 
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input less the losses over a given period of time. Therefore, the change in the average 
number of particles with respect to time is given by: 


dn 47rR 2 n. . v. . n 3n. . n 

tp J in J 1D J (3) 

dt 4/3 irR 3 t R/v. . t 

where t c is the confinement time. G. jt traM . The code is only a differencing solution of a 
one-dimensional boundary-value problem. Thus dn/dt can only be tied to the change in 
edge (boundary) density, as opposed to average density, so a parameter, F q = n/n e , was 
introduced to relate edge density to average density. Since this ratio is not fixed, F q is 
calculated in the code each time the dn/dt equation is solved (see below). A second 
parameter, the ratio of free-stream half-transit time with no potential well to the actual 
transit time, 


R/v. . 
F = inj 


(4) 


was introduced to allow the transit time to be factored out of the dn/dt equation. The 
final form of the dn/dt equation is then: 





Calculation of F 

n 


(5) 


The average number of particles between each pair of grid points is used to calculate the 
number of particles in that shell. This is summed over the entire grid and divided by the 
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total volume to arrive at n. A modification was necessary, however, for the electron 
calculation in the outer shell. In high density cases, the electron density change across the 
last grid space h = R - ^ = R(l— <r x >) can be very large (ca. 10 3 -10 4 x) and an algebraic 
average is not accurate. An exponential mean form was derived based on the assumption 
that the logarithmic gradient of the edge electron density varies inversely with the Debye 
length, A d , at the outer radius, R, d(ln n r )/dr = -nJX R , and with a parameter, /?, 
(typically P similar to 0.7), that characterizes the drop in E between the outer radius and 
the first grid point, cr^. When the electron density begins to increase rapidly in this 
region, EKXL will switch from an algebraic average to the following: 


2 ^ X 

n = n 

R (!-<!,>) R 


( 6 ) 


Electron and Ion Recirculation 


Although analytic forms have been derived for ion loss due to core upscattering 6 - i and 
from fusion reactions 8 , it has been found convenient to choose a constant ion recirculation, 
G., as an independent input to the code. Variation in G. affects only the ion injection 
current, I., as the circulating ion current, I c = LG., is fixed by the parameters of a given 
case. The code has the ability to adjust I. as needed (see below). 

The electron recirculation, G^, flUy be input at a constant value or it may be calculated in 
the code using the physics of electron confinement under mirror reflection (MR) and wiffle 
ball (WB) modes. This model is discussed in detail in other EMC2 technical 
notes/reports 9 ' 10 - 11 and its pertinent algorithms and defining equations are given in Figure 
6. The magnitude of B is a separate input parameter for the mirror/ wiffle model. The 
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auxiliary potential, introduced in KXL to model mirror reflection of the electrons, 
$ B = $mv 2 , in the cusp region, is input as the magnitude of pVB. 

Ion Losses Due to Fusion Reactions 


Fusion reactions of ions in the core and surrounding region will consume ions. This source 
of loss is important only at conditions of high core density and large ion energy. The total 
fusion rate is calculated in the code as part of the fusion and gain calculations (see below). 
This allows for easy inclusion of the ion losses due to fusion in the ion dn/dt equation. If 
the total fusion rate is DD fug fusions/sec, then the number of ions lost to fusion is 
2DD fus ions/sec. This loss term would go into the original dN/dt equation and the 
resulting dn g /dt equation would be: 


d !ll= 

3n . n j 

fn DD, 1 

e + fus 

1 

1 

dt ' 

F F 

n t 

l°i *.v 


t. 

tr&ns 


where F DD = 4/3 *R 3 A tranB is introduced to allow the transit time to be factored out. 


Ion Losses Due to Core Upscattering 


Ion losses due to core collisional upscattering can be accounted for by adding a loss term to 
the dn/dt equation for ions. Basqitfm the work of Rosenberg and Krall 6 , as modified by 
Bussard 7 , the upscattering time for ions, t up is calculated each time the dn/dt equation is 
solved, from 


1 


t 


up 


Z 4 1 

1.892xl0 13 f 2 (l+f) 2 A 1/2 R 


n(r) 

iw 575 * 


( 8 ) 


A factoring parameter, F u = t up /t tran#> was introduced into the dn/dt equation so that the 
ion losses due to upscattering are given by n e /(F u t trang )- This form is easily incorporated 
into the ion dn/dt equation which becomes 


dn 

e 

3n in- 

' 1 

1 ' 

DD . 

f ua 

1 

dt 

F n F t n ' 

G i 

F 

u _ 

F F 
r n DD 

^trans 


The fractional velocity spread, f = tfv/v, is a specified input parameter. If f = 0, no ion 
upscattering is calculated. This parameter can be estimated by comparison of scattering 
collision rates in the core, edge, and mantle regions 12 - 13 - h of the system, comparing 
transverse momentum vs. radial momentum exchange processes. Such estimates are not 
yet complete; when available the analytic forms that describe them can be added to the 
dn/dt equation, as above. 

Average Transit Time Calculations 

The data in the code is quite sufficient to calculate the ion transit time directly. Using the 
potential at each grid point and the input parameters E^, dE and dE perp , the average 
transit time across each grid point can be calculated and summed across the entire grid to 
arrive at the average ion transit time through the system. 

The calculation of electron tranril time, however, presents an interesting problem. The 

direct method used for ion transit time does not lend itself to an accurate calculation of the 

average electron transit time as the individual times vary greatly with the energy spread of 

the electrons. Electrons with a maximum dE will be turned early and have short 

perp J 

transit times while the electrons with maximum energy, E. .+dE, and no dE will pass 

OJ ’ inj ’ perp r 

through the bottom of the well and reach the core yielding much longer transit times. 
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Due to the dependence of F n on transit time, a small deviation in average transit time 
could continue to feed back into a large error in dn/dt. In short, the detailed 
computational capability of the code would lead to amplification of numerical errors in 
original data. In this circumstance, the internal data calculation was abandoned in favor of 
an analytic form based upon system parameters, which yields a more accurate and 
consistent averaged result and utilizes less computational time. As derived and discussed 
in another EMC2 note 14 (see for definitions of terms) the electron transit time is: 


t 


tot 


2R 

v . 
i n 



f<'w> 2 -<V 2 l 



<r 0 >+ 

to 

A 

V 

2 

f 1 

77' + 

1/2 

Y( 1_a q) 

+ f I 1 / 2 

1+ 5< 1 - a ,)j 

A 

V 


-1 


( 10 ) 


where a = <p /<p . , n' = 1-a (1 -ri ), rt = E /E , 

q *o' ^mm ' q' V’ 'o ao' o’ 

f 9 2 

<r w > = <r Q > + ^ <r Q > rj' and <r Q > is the core radius. 

This is a calculation of the maximum transit time. But the electrons occupy a spread in 
total and transverse energy. Across this energy distribution the individual transit times 
can vary markedly. Analysis of this variation suggests an exponential distribution of 
transit times across the band. Integrated averaging over this variation yields a correction 
factor to give average transit times of the average electron in terms of the total maximum 
transit time of equation (10). factor is found to be: 


1 n 


ft * 

trans 

2K7v“ 

in 




( 11 ) 
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Fusion Rates and Power Gain 


The KXL code contained a subroutine to calculate the approximate neutron rate (1/2 the 
fusion rate) for DD fusion. This subroutine was modified in EKXL to include fusion rate 
calculations for DT, D 3 He and p n B based upon the deuteron density (as the code does not 
currently accommodate different electron and ion injection radii, necessary for 
center-of-momentum, center-of-geometry placement of these ions). Calculation of 
system gross gain for all four fuel combinations was also added to the EKXL code, -'here 
gross gain is defined as 


b..E, 

1J fus 


G =- 


P. . 

mj 


n?(r) <r(E)v.(r)47rr 2 dr 


+ P: 


(12) 


The densities used for each fuel combination were: 

DT For DT fusion rates, the existing deuteron density, D, is divided in half and 
treated as 50:50 D:T 


D 3 He For a 50:50 D: 3 He distribution the total D + 3 He density will be 2/3 the 
existing D density for the same potential distribution. 

p n B For 50:50 p: U B t^otal density of p + U B is 1/3 the existing D density. 


Bremmstrahlung Losses 

Bremmstrahlung will be an important energy loss mechanism from the system at high 
density conditions. This was estimated from standard formulae 18 for bremmstrahlung 
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radiation density, integrated over the system volume to give total bremmstrahlung output, 
ignoring self-absorption by the plasma and any internal reflection from the boundary 
walls. This is discussed in a separate EMC2 technical note 17 . The resulting expression is 
used in the EKXL code to calculate bremmstrahlung gross output as 


P 


bremm 


1.69xl0 — 32 


n.(r)n (r)Z 2 VE (r)47rr 2 dr 

0 i e e 


(13) 


This is added to the injection power to determine net power input and system net gain. 


P. . +P. 

mj bremm 


(14) 


Ion Current Adjustment 


The problem with the anode control, discussed previously - allowing too many ions to 
build up in the system and thus requiring the maximum ion current to be within range of 
the final ion current - was addressed in EKXL with a different form of anode control that 
calculates the ion current density based on the core densities, the electron current density 
and G. „ The difference in core densities is given by 18 


n. - n = 3.32 x 10 


|8 K-Q 


X 


(15) 


where <p o is the core potential and is the potential at the bottom of the well (both in 
eV), r 0 is the core radius in meters and density is in 1/m 3 . The difference in densities at 
the target anode potential, is then 
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(16) 


nf -n' = 3.32 X 10° 




The change in density from n. to nj is then 


An. = n? - n. = 3.32 x 10 8 If™ 1 ? - + (An + n ) - n. 

ill 2 v e t' i 

r o 


An. = 3.32 x 10 8 - ** - (n. - nj + An e 


(17) 


Substituting equation (15) into equation (17) gives 


a ~*P ) 

An. = 3.32 x 10 8 - * - + An 

i l e 


(18) 


and 


3.32 x 10 8 G. 

Ji = ^T ( ^ ) + J 'a 


(19) 


where J. e are the circulating current densities. 


Automatic Operating Mode 

In the XL code a case achieved a steady-state solution or none at all. The KXL code has 
the ability to stop at any "snapshot" point by specifying the time duration in the input file. 
This allows a succeeding run to be started at the end of a previous calculation, using a 
special data file created when that calculation was completed, and to continue for the time 
duration specified in the new input file. This allows for the adjustment of time steps and 
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ion current necessary to achieve the final steady-state solutions previously unattainable 
with the XL code. The drawback is that the manual adjustments of duration time and ion 
currents require continuing operator input. In the EKXL code an "automatic" mode was 
introduced that eliminates the need for most manual restarts. In automatic mode, the code 
will run to the end of the specified time duration or a steady-state solution, whichever 
comes first, using a series of intermediate time steps (simulating the manual restarts 
required in KXL) automatically adjusting the length of each intermediate time step based 
upon whether or not the previous time step was successful. This is facilitated by the ’on 
current control and allows the code to run unattended as far as is possible without 
requiring operator intervention. 

IV. SAMPLE CODE RUN FOR EKXL V4.1 

Figure 7 shows a sample input file for EKXL version 4.1 as it is described above. This 
calculation ran in automatic mode with fixed ion current for a 15 keV well in a 92 cm 
device. The name of the log file, a sample of which is shown in Figures 8 a-b, and the 
dataset file, a sample of which is shown in Figure 9, are taken from the name of the input 
file, with the extension indicating the nature of the data. For example, the input file in 
Figure 7 is named BREM2b.INP. The log file in Figure 8 is named BREM2b.LOG and the 
dataset file in Figure 9 is named BREM2b.OUT. The last character of the file name is 
used to indicate the number of the run. In this example, the "b" indicates that it is the 
second run and the "Initializatfffl file" is set to "YES" in the input file. This calculation 
started where BREM2a.INP ended, and this is verified in the beginning of the log file 
(Figure 8 a) where it states that initialization file BREM2.TMP was read. 

The log file prints run diagnostics of the numerical calculation 4 to help the user in 
determining if and when problems occurred in the calculation. Figure 8a shows the 
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beginning of the log file, while Figure 8b shows the end. As seen in the figures, the log file 
also contains the transit times (Tau), the electron confinement time, the ion upscattering 
time (if this calculation is being used), the radius of the wiffie ball, <r b >, the actual ion 
and electron currents being used in the code, the fusion rates for DT, D 3 He and p n B, the 
plasma gain for all four fuels and the bremmstrahlung power (in watts) as calculated from 
the DD density. The end of the file indicates that the snapshot was saved and lists the files 
to which output was written. 

Figure 9 shows only the first portion of an output dataset file. This file contains the 
potential and ion and electron densities at each grid point for every successful snapshot in 
the run. The output is formatted and comma-separated which allows for easy import into 
spreadsheet or graphics software. A sample plot of this output data is shown in Figure 10. 
The beam profiles from the input file have been printed across the top of this chart and the 
fusion rates have been imported from the log file. 

Hundreds of cases have been run with the EKXL code modelling many areas of interest 
from basic experiments, to current SCIF parameters, to reactor regimes. Results of these 
code runs have been used extensively to better understand the physics involved in 
Polywell tm devices and also to improve the modelling, especially as it is applied in the 
code. A compilation of results of specific interest to the SCIF experiment is presented in 
some detail in another EMC2 technical report on these computer simulations 19 . The 
potential of this code, however, te not been fully realized and there still remain many 
areas of interest that have not been explored. 
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Figure 1. Maximum operational levels for successful Vlasov— Poisson "1— D" code runs. 

Figure 2. a. KXL case with R weighting factor. 

b. Same case with n weighting. Note improved core resolution and density. 


Figure 3. 


Core density vs. electron current for a fixed electron confinement. Slope is linear for 
currents greater 1000 A. 


Figure 4. 


Core density vs. electron current for varying anode heights. Slope approaches linear as 
anode approaches zero. 


Figure 5. Erroneous electron transit times as calculated in KXL. For all currents, electron energy 

is 100 keV. 


Figure 6. 
Figure 7. 


Figure 8. 


Figure 9. 


Figure 10. 


Sample input file for an EKXL v4.1 code run. 


a. Beginning of the log file created by EKXL for the run in Figure 7. 

b. End of the log file created by EKXL for the run in Figure 7. 


Beginning of the output dataset file generated by the EKXL run in Figure 7. 


Sample plot of the data in the EKXL output file, as shown in Figure 9. This graph 
was generated in SuperCalc v5.0. 



REFERENCES 


I Ensley, D.L., "Physic* Modeling Project", internal memorandum to Pacific— Sierra Research 
Corporation, October, 1986 

3 Ensley, D.L., "1ACCEL2/EACCEL2 User Guide", internal memorandum to Pacific— Sierra Research 
Corporation, April 1987 

3 McDonald, J., B. Goplen, "ACCEL User’s Manual", Mission Research Corporation Report 
MRC/WDC— R— 180, November 1988 

4 Smithe, D., "XL Reference Manual", Mission Research Corporation Report MRC/WDC— R— 226, 

January 1990 

5 King, K.E., R.W. Bussard, "Examination of Electrostatic Well Formation for Spherically Convergent 
Electron and Deuteron Beams", Energy/Matter Conversion Corporation Report EMC2— 0990— 02, 

September 1990 as presented at the American Physical Society Thirty— Second Annual Meeting of the 
Division of Plasma Physics, Cincinnati, OH, November 1990 

6 Rosenberg, M.(EMC2), N.A. Krall, "Ion Loss by Collisional Upscattering", Krall Associates Report 
KA— 90— 39, October 1990 

7 Bussard, R.W., "Core Collisional Ion Upscattering and Loss Time", Energy/Matter Conversion 
Corporation Report EMC2— 1090— 03, October 1990 

8 Bussard, R.W., "Fusion Lifetime Limits on Ion Upscattering", Energy/Matter Conversion Corporation 
Report EMC2— 0291-04, February 1991 

9 Bussard, R.W., K.E. King, "Electron Recirculation in Electrostatic Multicusp Systems: I — 

Confinement and Losses in Simple Power Law Wells", Energy/Matter Conversion Corporation Report 
EMC2— 0491— 03, April 1991 

10 Bussard, R.W., K.E. King, "Electron Recirculation in Electrostatic Multicusp Systems: II — System 
Performance Scaling of One— Dimensional ’Rollover' Wells", Energy/Matter Conversion Corporation 
Report EMC2— 0791— 04, July 1991 

II Bussard, R.W., K.E. King, "Phenomenological Modelling of Polywell tm /SCIF Multi— Cusp 
Inertial— Electrostatic Confinement Systems", Energy/Matter Conversion Corporation Report 
EMC2— 1191— 02, November 1991 as presented at the American Physical Society Thirty— Third Annual 
Meeting of the Division of Plasma Physics, Tampa, FL, November 1991 

12 Rosenberg, M.(EMC2), N.A. Krall, "Ion Loss by Collisions Outside the Core", Krall Associates Report 
KA— 91— 02, January 1991 

13 Krall, N.A., "Anisotropy Profiles for the Ion and Electron Plasma in the Polywell tm /SCIF", Krall 
Associates Report KA— 91— 09, April 1991 

14 Krall, N.A., "Core and Mantle Properties of Weibel Instability in the Polywell tm /SCIF", Krall 
Associates Report KA— 91— 10, May 1991 

15 Bussard, R.W., K.E. King, "Electron Transit Time in Central Virtual Anode Wells", Energy/Matter 
Conversion Corporation Report EMC2— 0291— 03, February 1991 

1® Book, D.L., "NRL Plasma Formulary", 1987 revised, Naval Research Laboratory, Washington, D.C. 

17 Bussard, R.W., K.E. King, "Bremmstrahlung Radiation Losses in Polywell tm Systems", Energy/Matter 
Conversion Corporation Report EMC2— 0891— 04, August 1991 

18 Bussard, R.W., G.P. Jellison, G.E. McClellan, "Preliminary Research Studies of a New Method for 


17 


Control of Charged Particle Interactions", Pacific— Sierra Research Corp. Report PSR 1899, 30 November 
1988, Final Report under Contract DNA001— 87— C— 0052, Appendix A, Sect, on Negative Potential Well 
Formation 


19 Bussard, R.W., K.E. King, "Computer Modelling of SCIF Experiment Performance", Energy/Matter 
Conversion Corporation Report EMC2— 0691— 02, June 1991 




Maximum Operational Levels 
for Successful Vlasov-Poisson "l-D" Code Runs 






Core density vs. electron current for a fixed electron confinement. Slope is 
linear for currents greater 1000 A. 




Core density vs. electron current for varying anode heights. Slope 
approaches linear as anode approaches zero. 
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Figure 5. Erroneous electron transit time as calculated in KXL. For all currents, 
electron energy is 100 keV. 
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Figure 6. Formulas used in the mirror/wiffle model for Gj calculation. 





Figure 8. a 


Operating mode is AUTOMATIC 

Auxiliary potential option is on. 

Initialization file read: brem2.tmp 

iter * 1, dumax ■ 3.21E-07, step « 1.00E+00 

Electron confinement time: 6.389E-07 

Wiffle radius rb: 1.936E-01 

1) time = 1 . 000E-06, determ = 1.000E+20 Tau = 6.100E-08 3.019E-06 

2) time * 1.263E-06, determ = 1.000E+20 Tau * 6.100E-08 3.027E-06 

3) time = 1.395E-06, determ = 1.000E+20 Tau = 6.100E-08 3.019E-06 

4) time = 1.972E-06, determ = 1.000E+20 Tau = 6.100E-08 3.022E-06 

5) time = 2.053E-06, determ = 1.000E+20 Tau = 6.100E-08 3.021E-06 

6) time = 1.527E-06, determ = 1.000E+20 Tau = 6.100E-08 3.021E-06 

Figure 8. b 

iter = 1, dumax = 6.82E-07, step = 1.00E+00 

Electron confinement time: 6.996E-06 

Wiffle radius rb: 2.692E-01 

787) time = 1.636E-04, determ = 1.000E+20 Tau - 6.113E-08 2.911E-06 

788) time = 1.663E-04, determ = 1.000E+20 Tau = 6.113E-08 2.911E-06 

789) time = 1.676E-04, determ = 1.000E+20 Tau = 6.113E-08 2.911E-06 

790) time = 1.736E-04, determ = 1.000E+20 Tau = 6.113E-08 2.911E-06 

791) time = 1.744E-04, determ * 1.000E+20 Tau = 6.113E-08 2.911E-06 

792) time = 1.690E-04, determ = 1.000E+20 Tau = 6.113E-08 2.911E-06 

iter » 1, dumax = 4.24E-07, step = 1.00E+00 

Electron confinement time: 6.995E-06 

Wiffle radius rb: 2.692E-01 

t (upscatter ) = 0.000E+00 0.000E+00 
Gun current 1 is: 1.800E+04 

Gun current 2 is: 1.800E+04 

Gun current 1 is: 1.800E+04 

Gun current 2 is: 1.800E+04 

Fusion Rate: DT*1.313E+14 DHe3=1.951E+10 DHe3 n-1.643E+ll pBll«5 . 434E+02 
Gplasma: DD=6.624E-09 DT=1.417E-06 DHe3=9. 550E-10 pBll«2 . 900E-18 
Pbrem: 1.059E+01 

Gicc*= 9 . 130E-06 No/Nc - 9.|^E+01 (<No>/Nc)*2- 1.378E+03 

##/ SNAPSHOT SAVED ### 

Restart file written: brem2.tmp 

Print files written: brem2b.out 

brem2b.tim 


Figure 8. 


a. Beginning of the log file created by EKXL for the run in Figure 7. 

b. End of the log file created by EKXL for the run in Figure 7. 


Start dates 07/16/91 Start time: 17:29:20.81 

Dataset : brem2b 

Time: 1.74E-04 sec 

Number of Grid Points: 200 

Deut Beam Neutron Count rate (#/sec): DD=1 . 479E+12 , DN=0.000E+00 


R <m) 

Pot (V) 

DD (n/scm3) 

Lambda (m) 

Nb (#/m**3 ) 
Elec 

Deut 

0.000000, 

-1 . 453E+04, 

1.718E+17, 

0.000E+00, 

6 . 601E+20, 

6 . 601E+20 

0.000434, 

-1 . 453E+04, 

1 . 718E+17 , 

0.000E+00, 

6 . 601E+20, 

6 . 601E+20 

0.000868, 

-1.453E+04, 

1.718E+17, 

O.OOOE+OO, 

6. 601E+20, 

6 . 601E+20 

0.001302, 

-1 . 453E+04 , 

1.719E+17, 

0.000E+00, 

6. 601E+20, 

6 . 601E+20 

0.001737, 

-1 . 453E+04 , 

1 . 720E+17 , 

0.000E+00, 

6. 601E+20, 

6 . 601E+20 

0.002171, 

-1 . 453E+04, 

1 . 722E+17 , 

0.000E+00, 

6 . 601E+20, 

6 . 601E+20 

0.002605, 

-1.453E+04, 

1 . 724E+17 , 

O.OOOE+OO, 

6.601E+20, 

6 . 601E+20 

0.003039, 

-1 . 453E+04, 

1 . 727E+17 , 

O.OOOE+OO, 

6.601E+20, 

6 . 601E+20 

0.003473, 

-1 . 453E+04, 

1 . 731E+17 , 

O.OOOE+OO, 

6. 601E+20, 

6. 601E+20 

0.003907, 

-1 . 453E+04 , 

1 . 736E+17 , 

O.OOOE+OO, 

6. 601E+20, 

6. 601E+20 

0.004341, 

-1 . 453E+04 , 

1.744E+17, 

O.OOOE+OO, 

6. 601E+20, 

6 . 601E+20 

0.004776, 

-1 . 453E+04 , 

1 . 753E+17 , 

O.OOOE+OO, 

6 . 601E+20, 

6 . 601E+20 

0.005210, 

-1.453E+04, 

1 . 765E+17 , 

O.OOOE+OO, 

6. 601E+20, 

6 . 601E+20 

0.005644, 

“1 . 453E+04 , 

1 . 780E+17 , 

O.OOOE+OO, 

6. 601E+20, 

6 . 601E+20 

0.006078, 

-1 . 453E+04, 

1 . 799E+17 , 

O.OOOE+OO, 

6. 601E+20, 

6. 601E+20 

0.006512, 

-1 . 453E+04 , 

1 . 821E+17, 

O.OOOE+OO, 

6 . 601E+20 , 

6 . 601E+20 

0.006946, 

-1 . 453E+04 , 

1 . 845E+17 , 

O.OOOE+OO, 

6. 601E+20, 

6 . 601E+20 

0.007380, 

-1 . 453E+04, 

1 . 870E+17 , 

O.OOOE+OO, 

6. 601E+20, 

6 . 601E+20 

0.007815, 

-1 . 453E+04 , 

1 . 893E+17 , 

O.OOOE+OO, 

6. 601E+20, 

6 . 601E+20 

0.008249, 

-1 . 453E+04 , 

1 . 908E+17 , 

O.OOOE+OO, 

6 . 601E+20 , 

6 . 601E+20 

0.008683, 

-1 . 453E+04 , 

1 . 909E+17 , 

O.OOOE+OO, 

6 . 601E+20 , 

6 . 601E+20 

0.009122, 

-1 . 453E+04, 

1 . 885E+17 , 

O.OOOE+OO, 

6. 601E+20, 

6 . 601E+20 

0.009574, 

-1 . 453E+04 , 

1.820E+17, 

O.OOOE+OO, 

6. 601E+20, 

6 . 601E+20 

0.010045, 

-1.457E+04, 

1.313E+17, 

O.OOOE+OO, 

5 . 760E+20, 

5 . 760E+20 

0.010543, 

-1.464E+04, 

7 . 805E+16, 

O.OOOE+OO, 

4 . 414E+20, 

4 . 414E+20 

0.011077, 

-1.468E+04, 

5 . 492E+16, 

O.OOOE+OO, 

3 . 689E+20, 

3 . 689E+20 

0.011655, 

-1.472E+04, 

4.027E+16, 

O.OOOE+OO, 

3 . 150E+20, 

3 . 150E+20 


Figure 9. 


Beginning of the output dataset file generated by the EKXL run in Figure 7. 


1) Elec 18.0e03 14.5e03 0.50e03 4.35e03 3.00e03 

2) Deut 18.0e03 lO.OOeOO l.OOeOO 1.71e00 leOO 

DD=1.479E+12 

Fusion Rate: DT=1.313E+14 DHe3=1.951E+10 DHe3 n=1.643E+ll pBll=5 434E+02 
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figure 10 Sample plot of the data in the EKXL output file, as shown in Figure 9. This 
graph was generated in SuperCalc v5.0. 




